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The interior structure of the Sun can be studied with great accuracy using observations of its 
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oscillations, similar to seismology of the Earth. Precise agreement between helioseismologi- 
cal measurements and predictions of theoretical solar models 1 has been a triumph of modern 



> ■ astrophysics. However, a recent downward revision by 25-35% of the solar abundances of 

<n : 

light elements such as C, N, O and Ne 2 has broken this accordance: models adopting the new 
abundances incorrectly predict the depth of the convection zone, the depth profiles of sound 
speed and density, and the helium abundance. 1 ' 3 The discrepancies are far beyond the un- 
certainties in either the data or the model predictions. 4 Here we report on neon abundances 
relative to oxygen measured in a sample of nearby solar-like stars from their X-ray spectra. 

H ' 

They are all very similar and substantially larger than the recently revised solar value. The 
neon abundance in the Sun is quite poorly determined. If the Ne/O abundance in these stars 
is adopted for the Sun the models are brought back into agreement with helioseismology 
measurements. 5 ' 6 

The role of the Sun as a fundamental benchmark of stellar evolution theory, which itself un- 
derpins much of astrophysics, renders the "solar model problem" one of some importance. The 
schism between helioseismology and models with a revised composition has arisen because abun- 
dant elements such as C, N, O and Ne provide major contributions to the opacity of the solar 
interior, which in turn influences internal structure and the depth at which the interior becomes 
convective. Uncertainties in the calculated opacities themselves appear insufficient to bridge the 
gap, 5 ' 7 and the propriety of the recent abundance revisions has therefore been questioned. 
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The revised abundances are demanded by new analyses of the visible solar spectrum that take 
into account convection and associated velocity fields and temperature inhomogeneities through 3- 
D hydrodynamic modelling, and that relax the assumption of local thermodynamic equilibrium 
for computing atomic level populations. 8 The solar abundances of C, N and O can be measured 
accurately based on their absorption lines. However, Ne lacks detectable photospheric lines in cool 
stars like the Sun, and the Ne abundance is therefore much less certain. The solar Ne content is 
assessed based on observations of neon ions in nebular and hot star spectra, and on measurements 
of solar energetic particles. 2 ' 9 Measurements are generally made relative to a reference element 
such as O; the downward revision of the solar O abundance therefore required a commensurate 
lowering of the Ne abundance for consistency. 2 Expressed as the ratio of the number of Ne atoms 
in the gas to those of O, this abundance is A^ e j Ao = 0.15, which is very similar to values adopted 
in earlier studies. 10 ' 11 However, it has recently been pointed out that the solar model problem could 
be solved were the true solar Ne abundance to be at least a factor of 2.5 times higher than recently 
assessed. 4,5 

We are motivated by the solar model problem to investigate Ne abundances in nearby stars. 
While not detected in the optical spectra of cool stars, emission lines of highly ionised Ne are 
prominent features of their X-ray spectra. 12 The Ne/O abundance ratio can be derived directly from 
the ratio of observed fluxes of the hydrogen-like and helium-like ions of O and Ne (see Methods 
for details); we adopt a slightly refined version of the method applied to the analysis of earlier solar 
X-ray spectra. 13 In this way, we have obtained the Ne/O abundance ratios for a sample of 21 stars 
lying within lOOpc of the Sun that have been observed by the Chandra X-ray Observatory using 
the High Energy Transmission Grating Spectrometer. 14 A representative X-ray spectrum, that of 
the Ml V star AU Mic, is presented in Figure[T] The stars studied, together with observed Ne and O 
line fluxes and derived Ne/O abundance ratios are listed in Table 1 . The Ne/O abundance ratios are 
illustrated as a function of the ratio of logarithmic X-ray and bolometric luminosities, Lx / Lboi — a 
commonly used index of stellar coronal activity — in Figured Since our star sample contains more 
objects toward higher L x / L bol , we have added Ne/O ratios for two stars of somewhat lower activity 
level (Procyon, an F5 subgiant, and e Eri, a K2 dwarf) from the literature. 15 There is no trend in 
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the Ne/O abundance ratio with L x /L boh and the error- weighted mean ratio is A Ne /A Q = 0.41. 

Our derived Ne/O abundance ratio is 2.7 times higher than the currently recommended solar 
value 2 but is consistent with the abundance inferred from helioseismology. 56 Solving the solar 
model problem by raising the Ne abundance alone would require a minimum ratio A Ne /A = 0.52, 
or 3.44 times larger than recommended. 5 However, raising the C, N, O and Fe abundances upward 
within their estimated uncertainty range of ±12% (adjusting them all together is not unreasonable 
because the recent downward revisions are correlated) would require an Ne abundance higher by 
only a factor of 2.5 5 ' 6 — quite within our estimated range. 

Extensive review articles on the coronae of the Sun and stars document evidence that coro- 
nal chemical compositions often appear different to those thought to characterise the underlying 
photospheres. 9 ' 16-18 The differences appear to relate to element first ionisation potentials (FIPs): 
the solar corona appears enhanced in elements with low FIP (< 10 eV; e.g., Mg, Si Fe), 9 ' 16 whereas 
much more coronally active stars appear to have depletions in low FIP elements and perhaps en- 
richments of high FIP elements (< 10 eV; e.g., O, Ne, Ar). 16-18 At first sight, this chemical frac- 
tionation might render direct interpretation of coronal abundance ratios in terms of the composition 
of the underlying star problematic. However, that we are seeing the same Ne/O ratio in a wide va- 
riety of stars sampling a large range of different coronal activity levels indicates that there is no 
significant fractionation between O and Ne in disk-integrated light from stellar coronae. This leads 
us to conclude that the results represent the true Ne/O abundance ratios of these stars. This con- 
clusion is bolstered by findings of a constant Ne/O abundance ratio, in good agreement with our 
value, for a small sample of single and active binary stars observed by XMM-Newton. 18 ' 19 In view 
of the consistency of the Ne/O abundance ratio in nearby stars, it seems likely that the solar ratio 
should be similar. 

There are no recent full-disk integrated light Ne/O measurements for the Sun; existing studies 
are instead based on observations of particular regions and structures of the solar outer atmosphere. 
We provide a tabular summary and bibliography of some of the different Ne/O estimates since 1974 
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as Supplementary Information to this letter. While the preponderance of estimates appear consis- 
tent with current assessments, the underlying observations do not sample photospheric material. 
The Ne/O ratio is in fact observed to differ substantially between the different observations. In 
particular, the highest measured Ne/O ratios based on X-ray lines are 2-3 times the accepted solar 
value and are compatible with the abundance ratio we find for nearby stars. These higher Ne/O 
ratios were obtained for hotter active regions 13 20 that are likely to dominate the solar full-disk X- 
ray emission. These measurements are the most directly compatible with the ones presented here 
based on full-disk integrated light X-ray spectra of stars. 

Similarly "high" Ne/O ratios have also been seen in 7-ray observations of flares, 3 He-rich so- 
lar energetic particle events, and in the decay phase of long duration soft X-ray events. 2123 The 7- 
ray measurements probe material that is irradiated by downward-flowing accelerated particles. 2122 
The particle beams penetrate through to the chromosphere which is likely more representative of 
the underlying photospheric material than coronal regions. 

Recent sophisticated models of heliospheric pickup ion and anomalous cosmic ray popula- 
tions from the local interstellar medium are also inconsistent with the current solar Ne/O abundance 
ratio, but could be reconciled by raising this to A Ne /A ~ 0.53, 24 which is similar to our values 
for nearby stars. This Ne/O ratio is somewhat higher than the mean from studies of ionized nebulae 
in the Milky Way and other galaxies, but falls within the scatter of results at solar metallicity. 25 

The implication of our study, then, is that the higher of the observed solar Ne/O abundance ra- 
tios are the ones representative of the underlying solar composition. This scenario is in accordance 
with our observations of Ne/O in nearby stars and reconciles solar models with helioseismology. 

Methods 

We use the resonance lines of H-like O and of H-like and He-like Ne to estimate the Ne/O abun- 
dance ratio. In hot (10 6 -10 7 K) coronal plasma these lines are formed predominantly by radiative 
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de-excitation of levels excited by collisions with thermal electrons. The flux, Fji, from such a 
transition j — ► i in an ion of an element with abundance A can be written as 



where Gji(T) describes the line emissivity — the product of the relative population of the ion in 
question and the excitation rate of the transition as a function temperature, T. The kernel $(T) — 
the emission measure distribution — describes the excitation power of the plasma as a function of 
temperature, which is proportional to the mean of the square of the electron density, N e , and the 



If we can choose O and Ne lines whose Gji(T) functions have very similar temperature 
dependence, an abundance ratio by number, A Ne /A Q , can be derived simply from the ratio of their 
observed line fluxes, Fq and F^ e , since all the temperature-dependent terms in Equation[l]cancel: 



An early study of Ne/O ratios in solar active regions 13 used the ratio of Ne IX \s2p 1 P 1 — > Is 2 1 S 
to O VIII 2p 2 P 3 / 2 .i/2 — > ls 2 S'i/2- This ratio does, however, have some significant residual de- 
pendence on temperature. 13 Here we remove much of this temperature dependence by addition of 
Ne X 2p 2 P 3 / 2 ,i/2 — > ls 2 5'i/ 2 ; our combined Ne Gji(T) function is G Ne = G NeIX + 0.l5G NeX - 
The resulting ratio Go /G^e is illustrated as a function of temperature in Figure|3] We have verified 
the small residual temperature sensitivity evident in the lower panel of Figure|3]to be negligible for 
our analysis by integrating the products of Go and G^ e with empirically-derived emission mea- 
sure distributions, $(T), for different stars, 12 ' 26 and for functions $ oc T a , with 1 < a < 4: the 
integrated emissivity ratio from these tests was / $Go dT j J ^Gn & dT = 1.2 ± 0.1. We conclude 
that the line ratio method is robust and the higher Ne/O abundance ratio found here will not be 
significantly changed through performing full emission measure distribution modelling. 




(1) 



emitting volume V, $(T) = N*(T) 




(2) 



We measured Ne and O line fluxes (listed in Table 1) from Chandra HETG X-ray spectra 
obtained directly from the Chandra public data archive (http://cda.harvard.edu). Final listed fluxes 
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for Ne X include small reductions (< 15% for 17 out of 21 or our stars, and 25-37% for the 
remainder) to account for a weak blend of Fe XVII at 12. 12 A. The Fe XVII 12. 12 A contribution 
was estimated by scaling the observed strengths of unblended Fe XVII lines at 15.26, 16.77, 17.05 
and 17.09 A (the strong 15.01 A resonance line was omitted to avoid potential problems with its 
depletion through resonance scattering) by their theoretical line strengths relative to the 12.12 A 
line as predicted by the CHIANTI database. Minor blending in the wings of the Ne IX 13.447 A 
line was accounted for by fitting simultaneously with the neighbouring weaker lines, comprised 
of a Fe XIX-XXI blend at 13.424 A and Fe XIX 13.465 A, following a detailed study of these 
features in the Capella binary system. 27 Since these blend corrections are generally very small, 
the uncertainties in these procedures have negligible (< 10% ) influence on the derived Ne/O 
abundance ratios. 
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Figure 1: A Chandra Medium Energy Grating X-ray spectrum of the Ml V star AU Mic. These 
data are typical of the Chandra MEG spectra upon which this study is based. Emission line features 
superimposed on the bremsstrahlung continuum are formed predominantly by radiative decay of 
transitions excited by electron impact. Prominent lines are labelled, including the lines of O VIII 

2j9 2 P 3 /2,i/2 -> Is 2 S 1/2 (18.97 A), NelXl^P! -> Is 2 % (13.45 A) and Ne X 2p 2 P 3 /2,i/2 -> 
Is 2 Si/ 2 (12.13 A) used here to determine the Ne/O abundance ratio. 
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Table 1: Spectral line fluxes and derived Ne/O abundance ratios for the stars anal- 
ysed in this study. Line fluxes were measured from the Medium Energy Grating 
(MEG) component of Chandra HETG spectra by line profile fitting using the Package 
for Interactive Analysis of Line Emission (PINTofALE) software 28 (freely available from 
http:hea-www.harvard.edu/PINTofALE). The effective collecting area of the instru- 
ment was accounted for using standard Chandra calibration products and techniques (see 
http://cxc.harvard.edu/ciao/ for details). Ne/O abundance ratios were obtained assuming the O/Ne 
line emissivity ratio of J §G dT j J §G Ne dT = 1.2 ± 0.1, as described in Methods. Stated flux 
and abundance ratio uncertainties correspond to la limits. 
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Figure 2: Derived Ne/O abundance ratios by number, An s /Ao, vs. the coronal activity index 
Lx/Lboi- Error bars represent quadrature addition of la uncertainties of line flux measurement. 
Also shown using hollow symbols are literature values 15 for the stars Procyon (F5 IV) and e Eri 
(K2 V) observed using the Chandra Low Energy Transmission Grating Spectrometer (LETGS) 
to better represent the lower ranges of coronal activity. The error-weighted mean Ne/O abun- 
dance ratio is A^ e j Ao = 0.41, or 2.7 times the currently assessed value 2 which is illustrated by 
the dashed horizontal line. The recommended value from comprehensive earlier assessments in 
common usage 10 ' 11 are also illustrated. 
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Figure 3: The temperature-insensitive O/Ne emissivity ratio. The upper panel illustrates the emis- 
sivities G , of the O VIII 2p 2 P 3 /2,i/2 -» Is 2 S 1/2 line, and G Ne of the Ne IX ls2p 1 P 1 -> Is 2 % 
and Ne X 2p 2 P 3 / 2 ,i/2 — > Is 2 S , i/ 2 lines combined as Gnc = G^eix + 0-l5GN e x- The lower panel 
shows the logarithmic ratio Go jG^ e . Emissivities are based on electron excitation rates and ion 
populations 29 compiled in the CHIANTI database, 30 as implemented in PINTofALE. 28 
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